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ABSTRACT 

Recent radio VLBI observations of the ~parsec-scale nuclear region of the 
narrow line Seyfert 1 galaxy NGC 4051 hint toward the presence of outflowing 
plasma. From available literature we have collected high-quality, high-resolution 
broadband spectral energy distribution data of the nuclear region of NGC 4051 
spanning from radio through X-rays, to test whether the broadband SED can 
be explained within the framework of a relativistically outflowing jet model. We 
show that once the contribution from the host galaxy is taken into account, the 
broadband emission from the active galactic nucleus of NGC 4051 can be well 
described by the jet model. Contributions from dust and ongoing star-formation 
in the nuclear region tend to dominate the IR emission even at the highest res- 
olutions. In the framework of the jet model, the correlated high variability of 
the extreme ultraviolet and X-rays compared to other wavelengths suggests that 
the emission at these wavelengths is optically thin synchrotron originating in 
the particle acceleration site(s) in the jet very close (few r g =GM-B}i/c 2 ) to the 
central supermassive black hole of mass Mbh- Our conclusions support the hy- 
pothesis that narrow line Seyfert 1 galaxies (which NGC 4051 is a member of) 
harbor a "jetted" outflow with properties similar to what has already been seen 
in low-luminosity AGNs and stellar mass black holes in hard X-ray state. 

Subject headings: galaxies: active — galaxies: nuclei — galaxies: Seyfert - 
galaxies: individual (NGC 4051) 



1. Introduction 



While it is generally agreed that the enormous energy output of active galactic nuclei 
(AGN) as well as X-ray binaries (XRB) is due to accretion of matter onto a compact object 
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(GM/[c 2 R] ~ 1, see, e.g., iReea Il984j for a review), the fate of the accreted material is 



somewhat uncertain. Most likely, a large fraction of this accreted matter crosses the event 
horizon of the black hole and is forever lost from the observer's view. However, the presence 
of relativistic outflows, often dubbed "jets" , observed to emanate from the very core of the 
AGNs and XRBs, show that some of the accreted material is in fact ejected from the inner 
accretion flow. Recent simulations show that such jets are likely Poynting flux dominated, 
i.e., only a very small fraction of the accreted material is ejected while most of the energy in 
the jets is carried by the electromagnetic fields. Launching of such relativistic outflows, their 
composition, and their association with the accretion inflow are currently topics of intense 
research. Observations indicate that nuclear radio-loudness of AGNs (parametrized by the 
radio-to-op tical luminosity ratio) tends to increase with decreasing Eddington-scaled nuclear 
luminosity ( lHoll2002t ISikora et al.ll2007l ). A similar situation is seen in XRBs accreting at low 
mass accretion rates where the radio loudness (in this case parame trized by radio to s oft X- 
ray flux ratio) tends to increase with decreasing X-ray luminosity (IFender et al.l 120031 ) . The 
radio emission from a radio loud AGN is generally attributed to synchrotron emission from 
a population of nonthermal particles in the jets and the lobes. But the emission mechanism 
in radio quiet AGNs su ch as Seyferts remains mu ch less well understood. Recent VLBI 
observations reported by iGiroletti &: Panessal (120091 ) have detected radio emission from the 
nuclear region of four very faint Seyfert galaxies (NGC 4051, NGC 4388, NGC 4501, and 
NGC 5033), where high brightness temperature (Tb ~ 10 5 -10 7 K) and (in some cases) 
resolved structure in the core favor a weak jet/outflow origin for these sources. 



Among these four faint Seyferts studied by Giroletti fc Panessal (120091 ) . NGC 4051, 
one of the "original" Seyfert galaxies studied by Seyfertl ( 1943h due to its strong nucl ear 



emission lines, is a narrow line Seyfert 1 galaxy (NLS1; see lOsterbrock fc Poggd Il985l for 
definition of NLSls) that has been well studied in almost all wavelengths from radio to 
the X-rays. Located in the Ursa Major cluster of galaxies, NGC 4051 is a fairly nearby 
(z=0. 002336 ±0.000004) spiral galaxy. Based on the latest reverberation mapping estimates, 
the AGN in t he heart of NGC 40 51 boasts a supermassive black hole (SMBH) of (1.7 ± 
0.5) x 10 6 M ( jPenney et al.ll2009l ). As mentioned above, recent VLBI observations of NGC 
4051 at 1.6 and 5 G Hz sh ow a complex, sub-parsec-scale compact core and other structure 



( foiroletti fc Panessal lioojjL Si milar results have a lso been obtained from recent VLA/ EVLA 



observations (IKing et al.l 1201 ll ; I Jones et al.l 120 111 ). The high brightness temperature of the 
compact core of NGC 4051 (> 1.8 x 10 5 K) along with the observed radio structures lead to 
the conclusion that the radio emission maybe be nonthermal in origin, and likely originating 
from an outflow/jet. 

From Gemin i Near-infrared Integral Field Spectrograph (NIFS) data of NGC 4051, 



Riffel et al.l (120081 ) have found a strong Hll Br7 emission at 2.1661 /zm in the nuclear spec- 
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trum (see their Fig. 1). Based on H2A2.1218/mi/Br7 line ratio they concluded that some 
contribution to the strong Br7 emission is from shocks or energy injection due to interaction 
of the H 2 emitting gas with the jet. While the strong Br7 line may be associated with 
interaction with the jet, it is quite unlikely that the IR continuum emission is dominated 
by the j et. Earlier high t ime resolution monitoring at 2.2yum continuum emission of NGC 
4051 by [Hunt et al.l ( 119921 ) does not show any rapid IR variability (relative to X-rays) which 
might be expected if the dominant fraction of the observed IR emission is an extension of the 
optically-thin X-ray emission (e.g. as r ecently found for the stellar black hole binary system 
XTE J1550— 564 bv lRussell et al.ll2010f). Simultaneous optical, IR and X-ray observations of 
NGC 4051 obtained by lDone et al.l (119901 ) also show little optical/IR variability compared to 
the X-rays, suggesting different origins for the X-ray continuum and optical/IR continuum 
emission. 

Compact, steady jets are observed from XRBs when they are in hard X-ray state, i.e., 
their X-ray spectrum is dominated by hard power law emission and large r.m.s. variability. 
In contrast, such steady, compact jets are not seen when XRBs are in the soft X-ray state 
characterized by pre dominantly therma l spec tr um from the accretion d i sc and low r.m.s. 
variability. See, e.g., iHoman fc Bellonil (120051 ); iRemillard fc McClintockl (120061 ). for exten- 
sive definitions and review of X-ray states in XRBs, and Fender for a review of XRB 
jets. It is still not clear which physical mechanism(s) trigger transition between these states, 
as well as the associated change in disk-jet coupling. Such state transitions and appear- 
ance/disappearance of radio jets are now regularly seen in transient X-ray binary systems 
(where epochs of activity typically span ~month or less), thanks to advances in instrumenta- 
tion and coordinated, quasi-simultaneous multiwavelength observations over the past decade. 
Such changes however cannot be seen in AGN because the associated timescales can easily 
exceed the human lifespan, and advances can only be made via studying ensemble proper- 
ties of AGN. Indeed, recent works show that there exists a tight relationship connecting the 
radio luminosity, X-ray luminosity, and black hole mass. This relationship, often referred t o 
as the "fundamental plane" of black hole accretion (IMerloni et al.ll2003l ; iFalcke et al.l 120041 ). 
not only holds across the stellar mass black holes in hard state and low-luminosity AGNs 
(LLAGN), but also for sup ermassive black holes with direct, dynamical mass- measurements 
(i.e. the M — a AGNs; see lGultekin et al.l 120091 ). and extends to blazars as well (Plotkin et 
al. in prep.). 



McHardy et al.l (120041 ) have argued that NGC 4051 and the stellar mass black hole 



Cygnus X-l in soft state have very similar Fourier power density spectra, and therefore NGC 
4051 may be a soft state AGN accreting at ~30% L Edd . Based on the bending frequency 
of the Fourier power spectrum of NGC 4051, these authors had deduced a central black 
hole mass of 3±l x 1O 5 M , which was consistent with the best estimated reverberation 
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mass of 5^3 x 1O 5 M at that tim e. However, the mos t recent, revised reverberation mass 



estimate of 1.73±°;g 



x 10 M by iDenney et al.l ( 120091 ) is ~5.7 times larger, consequently 
lowering the bolometric luminosity to ~5% Lead- Moreover, as discussed in detail in their 
discussion section, the assumed X-ray to bolometric conversion factor of 27 is quite likely 
an overestimate. On the other hand, another commonly used empirical rel ationship to 
estimate the bolometric lu minosity of similar AGNs is Lboi ~ 9ALa (5100 A) ( Kaspi et al. 



2000l ; iPeterson et al.l 120041 ) . Given NGC 4051 's host galaxy starlight subtracted averag e 



nuclear flux density at 5100 A is (4.5±0.4) x 10~ 15 erg s" 1 cm -2 A" 1 (IDennev et all 120091 ) 



its bolometric luminosity is roughly « (1%)(D/10 Mpc) 2 (M BH /l-7 x lO b M ) _1 L Ed d- Both 
of the above estimates of the bolometric luminosity are closer to Galactic black hole binaries 
in hard state than in soft state. It is also worth noting that Cygnus X-l is act ually never in 



a "classical" soft state, b ut rather it is always somewhat transitional (see, e.g.. iNowak et al. 



20051 . 1 Wilms et all 120061 ). Thus it is unclear whether NGC 4051 is indeed a soft state AGN, 
or in hard state, or perhaps in a transitory state between canonical hard and soft states. Can 
XRB state definitions be straightforwardly carried over to AGNs? Such are the questions 
which motivate the current work. 

Since jets are known to emit over a broad range of the electromagnetic spectrum span- 
ning from radio through X-rays, it is plausible that at least some of the observed nuclear 
emission from NGC 4051 might also originate in the jet. In this work we have collected 
high quality, high-resolution nuclear SEDs of NGC 4051 from the literature, to study the 
broadband (radio-X-ray) SED and test whether the observed nuclear emission can origi- 
nate in a mild l y rela t ivistic jet out fl ow. W e use the accretion disk+jet model developed by 
Markoff et al.l (120051 ). iMaitra et al.l ( 12009al ) to model the broadband continuum. The model 
has been successfully applied to explaining broadband SEDs of several black hole X-ray bi 



naries in hard state, as well as that of the LLAGNs Sagittarius A*, M81* (see, e.g.. iMarkofl 



2010l . and references therein). These previous results suggest that despite a difference of ~5- 
6 orders of magnitude in mass, the accretion/ejection processes in stellar mass black holes 
(~1OM ) and supermassive black holes in the nearby low-luminosity AGNs (~1O 6 ~ 7 M ) 
appear quite similar when comparing sources with similar Eddington scaled luminosity. 

This is the first time we apply this jet model to a radio-quiet Seyfert galaxy to test 
whether the broadband emission originates in a collimated jet outflow. The sources from 
which we obtained the multi-wavelength SED data are described in §2j and a brief summary 
of the model parameters that best match the data is given in §31 Implications of the modeling 
results and conclusions drawn from the study are discussed in §H 
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2. The broadband SED of NGC 4051 

The nuclear region of NGC 4051 has been observed with many different ground-based 
as well as space-borne telescopes spanning from radio through X-rays, at the highest spa- 
tial resolution available in each of these wavelength ranges. The high angular resolution is 
of crucial importance for isolating the central, nuclear emission from that of the surround- 
ing galaxy. This is especially important in optical and IR wavelengths where contributions 
from non-AGN components such as the host galaxy's stellar population, ongoing star for- 
mation near the nuclear region, and dust can easily dominate the observed emission if the 
photometric aperture is not small. 



2.1. X-ray 



The X-ray flux from NGC 4051 is known to be highly variable on a timescale of day s 



(see, e.g., iMcHardv et all 12004 : IPonti et al.l 12009 : Breedt et al.ll2010l ; IVaughan et aUEOllt ). 
Like most NLSls, the 2-10 keV X-ray continuum SED of NGC 4051 is well modeled by a 
power law. A narrow iron line is seen near 6. 4 keV, and the pr e sence of curvature at higher 
energies hint towards presence of reflection ( IPonti et al.l l2006t iMiller et al.l 120101 ) . Similar 
to the X-ray spectra of many other NL Ss, there exists an excess of soft photons in addition 
to the power law at <2 keV (see, e.g., iBoller et al.lll996l . also see S2J. In this work we are 
interested in the power-law continuum only and take the power-law photon index (r) and 
normalization from a rec ent simultaneous r adio ( VLA/ EVLA) -X-ray [Chandra) campaign of 
NGC 4051 published by iKing et al.l ( 1201 ll ). We chose this set of X-ray observations because 
they were performed simultaneously with the radio observations at 8.4 GHz. The solid black 
lines in Fig. [T] show the power-law continuum measured from this Chandra data set. As 
discussed in §4J the hard X-ray (>10 keV) variability may be smaller than the soft (2-10 
keV) X-ray variability. The refore we have also added two Suzaku/PIN spectra of NGC 4051 
from the data presented by IMiller et al.l (|2010|). While NG C 4051 has not been detected in 
the Fermi/LAT data so far (0.1-100 GeV; Abdo et al. 201oL it has been detected in the time 
averag ed Swift/BAT (14-195 keV) data. The BAT spectrunfl is also included in Fig. [TJ Fits 
to the high-resolution Chandra HETGS spectra of NGC 4051 (King et al., in prep.) suggest 
that any absorption intrinsic to the source is negligible, therefore to correct for absorption 



we us e a column density of Ah = 1.3 x 10 cm in the direction of NGC 4051 ( lElvis et al. 
1989h . 



1 Obtained from |http://heasarc.gsfc.nasa.gov/docs/swift/results/bs9mon/bsind.php?number=85 
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2.2. UV 



The Extreme Ultraviolet Explorer (EWE) data are from lUttley et al.l (l2000h. and the 



corre s ponding; error b ar rep resents the EUV variability range (also see Fig. 7 of iKaspi et al. 



2004 ). lUttley et al.l (l2000f ) noted that the EUV variability is comparable to the X-ray vari- 
ability, and that the EUV and X-rays were simultaneous to within 1 ksec, which suggests a 
common origin for both EUV and soft X-ray emission. 

The far-ultraviolet data in the wavelength range of 900-1180 A ar e taken from t he Fa r 
Ultraviolet Spectroscopic Explorer (FUSE) observations presented by IKaspi et al.l ( 12004 ). 
Even thoug h the three FUSE observations were separated by as much as ~ one year, the 
analysis by IKaspi et al.l ( 12004) suggests little or no (<10%) variation in the FUSE spectra. 
As discussed above, we are interested in the continuum emission only, and therefore use 
the average flux-density measurement s in wavelength ran ges that are free from emission and 
absorption as presented in Table 2 of IKaspi et al.l ( 12004 ). 



We have also included HST/STIS data of NGC 4051 published bv lCollinge et al.l (|200l|). 
The binned STIS spectrum presented in their work shows that the continuum flux is roughly 
constant at ~ 2 x 10~ 14 erg cm -2 s _1 A -1 over the observed wavelength range of 1200-1700 
A. These HST observations were made simultaneously with Chandra, and e ven though the 



Chan dra light curve showed large variations in the X-ray flux (see, e.g, Fig. 2 of lCollinge et al. 



20011 ) . no variability could be detected in the STIS data. The lack of UV variability strongly 
contrasts with the high X-ray variability and suggests that the UV and X-rays are most 
likely coming from different emission mechanisms. 



2.3. Optical 

Contribution of non-AGN components, e.g. light from stars as well as dust, starts 
becoming important in optical and IR, and extreme caution must be taken not only to use 
highest resolution available to observe the nuclear region, but also to model and subtract the 
non-AGN contributions within the aperture. Readily available optical and IR fluxes listed, 
e.g. in the NED, are therefore not usable for this work since the apertures listed in NED are 
typically large and there is no easy way to estim ate the non-AGN contribution. Instead we 



use the results obtained by iDenney et al.l ( 120091 ). who presented new, optical spectroscopic 



and photometric observations of NGC 4051 from a densely sampled reverberation mapping 
campaign, using data from five different observatories (1.3m McGraw-Hill telescope at MDM 
Observatory, 0.7 and 2.6m telescopes of the Crimean Astrophysical Observatory (CrAO), 
2m MAGNUM telescope at the Haleakala Observatories, and the 0.4m telescope of the 
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University of Nebraska(UNebr)). They then subtr act the host g alaxy starlight contribution 
to the continuum flux following the prescription by lBentz et al.l (120091 ) . The estimated AGN 
continuum flux at 5100 A for NGC 4051 is (4.5 ± 0.4) x 10" 15 erg s" 1 cm" 2 A" 1 after 
subtracting host galaxy starlight contribution. The host galaxy starlight contribution is 
about twice that of the AGN. However the optical flux is much le ss variable than the X- 
ray and EUV. The excess variance in the 5100 A measurements of iDenney et al.l (120091 ) is 
0.09. A recen t anal ysis of last twelve years of X-ray and optical variability in NGC 4051 by 
Breedt et al.l (120101 ) also concluded that the variability of the AGN in optical wavelengths 
(where the fractional variability, F var , is ~12-26%) is much smaller than the the variability 
in X-rays (F var ~ 49%). 

2.4. Infrared 



For I R we used the Sp itzer mid-infrared (MIR; 5-15 /zm) spectrum of NGC 4051 pub- 
lished by ISani et al.l (120101). an d also adaptive optics Integral Field Unit spectroscopic data 
reported by iRiffel et al.l (120081 ) which was obtained with the NIFS mounted on the Gemini 
North telescope. Three representative continuum fluxes at 2.1, 2.2, and 2.25/im from the 
nuclear (central 0.1"x0.1") region obtained by the NIFS are shown by the cyan points in 

Fig.m 



The Spitzer spectrum (also see Fig. 2 of ISani et al.l 120101 ) shows strong polycyclic aro- 
matic hyd r ocarb on (PAH) features due to ongoing star formation activity. According to 
Sani et al.l (120101 ). the MIR continuum may also be primarily due to thermal emission by 
AGN- heated dust. Therefore contribution of direct AGN emission in MIR appears to be 
small. The slit width of the IR spectrograph is 3.6 arcsec, which at the distance of NGC 
4051 implies a linear scale of ~175 (D/10 Mpc) parsec. The Spitzer/IKS slit width is a fac- 
tor of 18 larger than that of the HST/STIS slit width (0.2 arcsec) and the former therefore 
includes significantly larger circumnuclear regions. The Spitzer data therefore do not put 
strong constraints on our broadband modeling. However they put upper limits on the direct 
(accretion/ejection origin) AGN emission in the MIR wavelengths. The Gemini NIFS nu- 
clear spectrum with an angular coverage of 0.1" x 0.1", is expected to contain a significantly 
larger fraction of AGN-to-non-AGN (starburst, dust) emission compared to the IRS. 

2.5. Radio 



The 8.4 GHz radio data are taken from iKing et al.l (120111 ). where the authors present 
simultaneous radio and X-ray monitoring observations of NGC 4051. During this campaign 
NGC 4051 was observed 6 times between 2008 December 31 and 2009 July 31 simultaneously 
us ing both VLA/ EV LA GHz and Chandra. Details of the radio data analysis are presented 



in 



King et all (1201 If ). 
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The UV, optical and IR da ta were dereddened assuming E(B-V) =0.013 magnitudes 
and using the extinction law of ICardelli et al.l (119891 ) . We note that no cosmological K- 
correction was made to the observed continuum fluxes. Given the small redshift of NGC 
4051 (z=0.0023), any K-correction is negligible and uncertainties in absolute flux calibrations 
are much larger than this correction. Table [1] gives the wavelength range and observation 
dates as well as references to the published works we have used to construct the broadband 
SED. 

3. Modeling the SED 

The continuum j et emission in the out flow- dominated model we have used here is based 
on previous works by iMaitra et al.l (j2009al ) and iMarkoff et al.l (120051 ) . We refer the reader to 
these works for the details of the model. The main parameters that determine the properties 
of the jet in this model are the input jet power (Nj), electron temperature of the relativistic 
thermal plasma entering at the jet base (T e ), the ratio of magnetic to particle energy density 
(a.k.a. the equipartition factor k), physical dimensions of the jet base (assumed to be 
cylindrical with radius r and height h ), and the location of the point on the jet (z acc ) 
beyond which a fraction of the leptons are accelerated to a power-law energy distribution 
(dN/dE ~ E~ p ). Nj, parametrized in terms of the Eddington luminosity, determines the 
power initially input into the particles and magnetic field at the base of the jets. Beyond the 
jet-base, we assume that the jet plasma (a) accelerates mildly due to longitudinal pressure 
gradient with its bulk Lorentz factor asymptotically reaching a value of 2-3 far downstream 
along the jet, (b) expands laterally with its initial sound speed, and (c) cools due to adiabatic 
expansion. Once the initial conditions are specified, then the bulk speed of the jet plasma, 
magnetic field, and pa rticle density are computed by solving the adiabatic, relativistic Euler 
equation (jFalckdll996l ). We note that in AGN, a combination of higher powers, confinement 
and perhaps stronger magnetic field pressures may give higher bulk Lorentz factors. 

The a ccretion disk is modeled as a multi-color blackbo dy with radial temperature profile 
TocR -0 - 75 (jShakura fc Sunyaevlll973l ; iMitsuda et al.lll984j ). and is parametrized by T in , r in , 
the temperature and radius at the inner edge of the disc. The thermal photons from the 
accretion disc are also included in the photon field of the jet for Compton scattering, although 
due to Doppler redshifting in the jet frame, this external Compton emission is usually much 
smaller than the synchrotron self-Compton. 

Given the various uncertainties involved, e.g., in our understanding of the cause of the 
large X-ray variability, non-simultaneity of data from different wavelengths (except for the 
VLA/EVLA data and the Chandra data, which were obtained simultaneously), contribution 
to the optical/UV flux from non-AGN components, dust heating, etc., a detailed statisti- 
cal comparison between the data and the model would yield errors which will significantly 
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underestimate the systematic uncertainties. In other words, the uncertaintites involved are 
much larger than simply the deviation of model prediction from the data. Therefore the 
best-match model here is 'x-by-eye', and is shown by the line labeled 'Total' in Fig. [TJ In 
Fig. [2] we show the model variability when one of the parameters was varied by ±10% and 
±20% from the best-match model while the remaining parameter values were fixed to those 
of the best-match model. 



We adopted a distance of 10 M pc (ISani et al. 120101) and the recently revised reverbera 



tion mass estimate of 1.7xl0 6 M ( Penney et al. 20091 ) for NGC 4051. The outer radius of 



the accretion disc is not constrained by the data and was arbitrarily set to a value of 10 8 r g . 
The remaining parameters were obtained from the best-match to the data. 

Assuming that the radio originates in an optically thick, stratified jet, a comparison of 
the radio and Spitzer/IKS fluxes immediately shows that the Spitzer fluxes are at least an 
order of magnitude higher than the expected jet flux. The large slit-width of the IRS (3.6") 
is primarily responsible for this. The high spatial-resolution Gemini/NIFS spectrum, shows 
that the AGN contribution is much less. 

The soft X-ray conti nuum, obtained w ith Chandra simultaneously with the 8.4 GHz 
VLA/ EVLA observations (King et al.ll201ll ). shows much more variability both in the power 
law index and normalization. Such variations could arise due to rapid changes in the parti- 
cle energisation/acceleration/cooling processes occurring in s mall "knots" in t he po st-shock 
region of the jet, e.g., as seen in the blazar BL Lacertae by IMarscher et al. J2008f ). A full 
understanding of particle acceleration in such jets is still lacking; here we use a particle 
energy distribution index p=2.7 which matches the observed >10 16 Hz SED quite well if 
the observed EUV and X-ray emission is largely due to optically thin synchrotron. The jet 
inclination was taken to be 50°. While not strongly constrained by the data, a much lower 
value (<40°) of the inclination would cause a larger fraction of the Doppler-boosted syn- 
chrotron photons from the outer parts of the jet to overestimate the radio flux. Inclinations 
greater th an ~60°are similarly not p referred. We note that this inclination is same as the 
results by IChristopoulou et al.l (119971 ). who modeled the narrow-line region of NGC 4051 as 
an outflowing biconical region of inclination 50°, thus suggesting that the relativistic jet as 
well as the much slower outflowing narrow-line region are co-aligned. Since we assume a 
symmetric model where the jets are perpendicular to the accretion disc, this means that the 
inclination to the disc normal is also 50°. 

Since the radio data do not show any hint of low-energy cutoff, the minimum length of 
the jet was set to 0.3 p arsec, which is the upper li mit for the source radius of the central 
component detected by iGiroletti fe Panessal (120091 ) using VLBI. Also, while searching for 
the best-match model more weight was given to the VLA A configuration data point of 
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King et al.l ( 1201 ll ) since that observation had the highest spatial resolution as well as the 



smallest error bar of all the radio observations presented in that paper. 

The equipartition factor (k) is not strongly constrained by the data, although we could 
not obtain good matches either for k < 0.5 or for k > 10, pointing toward a near-equipartition 
model. We therefore set k=l. In addition to the jet, an accretion disk with ri n =130 r g and 
Ti n =2.8xl0 K is required to model the SED. This however assumes that the entire optical 
and UV emission is from the accretion disc. As discussed in §H (and perhaps as evidenced 
from the PAH lines in the Spitzer/IKS spectrum), some of the UV emission could be from 
an actively star-forming nuclear star cluster. Some of the optical/UV radiation could also be 
due to reprocessing of the X- rays from the jet ba se by the outer accretion disk, as suggested 



by the time-lags reported in iBreedt et al.l ( 120101 ) . Therefore the inner disc radius obtained 



here is more realistically an upper limit on the actual radius. 

The pre-shock synchrotron is the emission from the base/launching region of the jet. 
While the data does not constrain it directly in this case (see, e.g., the solid orange curve 
in Fig. [TJ, the physical conditions at the base determine the nature of the flow farther 
downstream along the jet in the model. In order to produce the observed radio emission 
(which arises in the optically thick, post-shock synchrotron regime) and the X-rays (which 
in our model are originating mainly in a region close to the base where the a fraction of 
the thermal leptons are accelerated to a non-thermal powerlaw distribution), only certain 
physical conditions and geometries of the base are allowed, and these choices (r =4 r g , 
h /r =l, T c =4xl0 10 K, z acc =5 r g ) lead to the observed pre-shock synchrotron emission. An 
input jet power of 1.7xl0 -2 LEdd, initially input into the particles and magnetic field at the 
base of the jets, is needed to match the broadband SED. The inverse Compton component 
(solid blue line in Fig. [TJ) is also not strongly required by the data, however it is computed 
self-consistently by the model and we include it in the plot to show that the X-ray emission 
is indeed dominated by synchrotron radiation and not by inverse Comptonization. We also 
tried models where the high energy is primarily dominated by inverse Compton scattering. 
However, within the framework of jet models, an IC-dominated model requires very extreme 
and rather unphysical parameters (e.g., r <r g , and h /r >20) for the X-rays, and is still 
unable to model the EUVE flux. Therefore we favor a synchrotron origin for the EUVE and 
soft X-ray flux. 



Discussion and Conclusions 



It has long been known that the so-called radio-quiet S eyfert ga laxies are not ent irely 



radio silent (see, e.g.. lUlvestad fc Wilsonlll989UWilsonlll99ll . and also lHo fc Pengil200ll for a 
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discussion). Recent VLBI observati ons of some of the faintest radio cores in Seyfert galaxies 
have suggested the presence of jets (IGiroletti Panessall2009l ). In this work we have studied 
the broadband nuclear emission from NGC 4051, a narrow line Seyfert galaxy, where the 
presence of a jet is suggested not only from radio VLBI, but also from observations at other 
wavelengths, e.g. presence of s hocked Br7 emis sion seen in high-resolution Gemini/NIFS 
spectrum of the nuclear region ( iRiffel et al .1120081 ). 



The modeling presented in this work suggests that once contributions from the host 
galaxy are properly accounted for, the broadband nuclear SED of the narrow line Seyfert 
galaxy NGC 4051 can be reasonably well described by emission from an accretion disk plus 
a mildly relativistic jet outflow. A comparison of the best-match model parameter values 
obtained for NGC 4051 with XRBs in hard state and the LLAGNs Sgr A* and M81* is 
given in Table |5J The jet parameters for NGC 4051, viz. input power (Nj), geometry of 
the launching region (r , h ), the temperature of the thermal plasma (T e ), equipartition 
factor (k), particle (lepton) energy distribution index (p) are in the same ballpark regime 
for stellar mass black hole X-ray binaries in hard state, even though the masses are different 
by 5 orders of magnitude. However, the location of the region beyond which particles are 
accelerated (z acc ) appears to be closer to the central engine for NGC 4051 when compared to 
the the LLAGNs. Broadband SED of more NLSls would be needed to confirm whether NGC 
4051 is an anomalous source in this respect, or if this is a generic feature of NLSls. From 
our modeling it appears that the jet axis m ay be co- aligned with the outflowing narrow-line 



region (see, e.g. 



Christopoulou et al.lll997l ). 



Even though there is not enough high-quality high energy (>20 keV) data to constrain 
the efficiency of particle acceleration in the post-shock jet, the Swift/BAT data suggests 
that the high energy cutoff, if at all present, is at >100 keV. Similarly it is not possible 
to strongly constrain the equipartition factor in this case, but models near equipartition 
or slightly magnetically dominated models are preferred. The location of the region where 

5 r v in t his case) is still a free parameter in the model. Recent 



acceleration start s (z a( - n ~ 



calculations by iPolko et al.l (120101 ) have indicated that it is possible to self-consistently de- 
rive the location of this acce leration region with i n the self-similar magnetohydrodynamical 
approximation proposed by IVlahakis fc Konigll ( 120031 ). Once gravity is incorporated into 
this new calculation, we will be able to reduce the number of free parameters significantly 
in future versions of the jet model (Polko et al., in prep.). 
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4.1. Variability 



The high variability of th e X-rays and the extreme- UV, and the fact that they are 
simultaneous to within 1 ksec (jUttley et al.ll2000l ) suggests a common origin for these two 
bandpasses. The variability of the optical and lower frequencies is much smaller. In the 
framework of our model, the EUV and X-rays originate in a compact jet close to the central 
black hole. The observed high variability is expected both because of the small size of the 
emission region and relativistic beaming due to bulk motion of the jet plasma. The optical, 
on the other hand, originates from a much larger region in the accretion disc, and therefore 
exhibits much lower variability. 

The energy distribution index of the shocked leptons (p = 2.7) o btained here is somewha t 
steeper than the values predicted e.g. by diffusive shock process (IHeavens fc Drurylll988l ). 
It is however possible that for NGC 4051 the EUVE and X-ray bandpasses are at a higher 
energy than the cooling break energy, where the particle energy distribution steepens from 
p to p + 1 due to enhanced cooling. If the cooling break is indeed at a frequency lower than 
the EUVE, an immediate effect of this would be to move the optically-thick-to-optically- 
thin transition, currently in the ~10-15 micr on range, to submm wavelengths (where this 
break has been seen for the LL AGNS Sgr A 



Falcke fc Markoff 2000; Melia k Falcke 2001 



Maitra et alJl2009b| . and M81* jMarkoff et al.ll2008| ). Since the location of the cooling break 
is dependent on time as well as particle injection/acceleration rate, some of the observed 
variability could also be attributed to cooling. 



Ponti et al.l (120061 ); iMiller et al.l ( 120101 ) have noted that above >5 keV there may be a 
tendency for the rms variability to decrease with increasing energy. X-ray flux in our model 
originates in the optically thin synchrotron from a distribution of shocked particles beyond 
z acc . The spatial extent of such shocks may be small, and the X-ray variability could be due 
to rapid changes in particle energisation, e.g, perhaps due to magnetic reconnection events. 
While the soft X-rays (< 10-20 keV) originate mainly due to this "post-shock" synchrotron, 
it is possible that inverse Compton scattering (of synchrotron photons produced at the base 
of the jet (SSC) as well as that of the soft photons from the accretion disc) becomes dominant 
at >20 keV. Assuming a reasonably constant mass accretion rate by the disc, and a constant 
fraction of mass channeled into the jets, the variability of this IC component would be small, 
and would largely reflect the (more slowly occurring) changes in mass accretion rate. 



The simultaneous Chandra and VLA data presented by King et al.l (120111 ) suggest that 
during their 2008-2009 campa ign NGC 4051 sho ws an anti-correlation between X-rays and 
radio. Another recent work by I Jones et al.l (120111 ). using RXTE and VLA observations made 
during 2000-2001, found a weak positive correlation betwee n radio and X-ray fluxes. The 
origin of the anti-correlation observed by King et al.l (120111 ) is still not clear; however this 
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may be due to the following two reasons: (a) given the uncertainties in the bulk speed of 
the jet plasma, the simultaneous (as observed by us) radio and X-ray fluxes may not be 
causally connected, and (b) as multiwavelen gth monitoring results of Sgr A* show (see, e.g. 



Eckart et al.ll2006l ; lYusef-Zadeh et al.l 120081 ). there seems to be no tight correlation between 



the X-ray and radio flares. In Sgr A* the optically thin NIR/X-ray flares most likely trace 
particle acceleration/(re)energization and cooli ng very close to the black hole, which only 
marginally affects the optically thick radio flux (IMarkoff et al.l 120011 ). Changes in radio flux 
on the other hand ar e mainly driven by adiabatic expansion of overdensities in the outflow 
( IMaitra et al.ll2009bl ). It is also possible that NGC 4051 is in a transitional phase between 
hard and soft states, where the radio emission is in the process of getting quenched. While 
comparing with XRBs one must however keep in mind that the observed "quenching" of 
XRB jets may not mean that the jets have actually turned off. It is possible that XRBs 
jets may only appear quenched because they have become quite weak and the sensitivity is 
not enough to detect them. The presence of a separate mode for the disk-jet coupling (i.e., 
separate from the cou pling which gives rise to the fundamental plane relation) operating in 
NGC 4051 during the 



King et al.l (120 111 ) observations cannot be ruled out either. 



4.2. The inner-edge radius 



Our modeling suggests that the inner edge of the accretion disk is ~130r g . Given that 

10~ 2 LEdd, a large r it1 is not unreasonable. The profil e of the 



the model luminosity is ~iu " L/Edd; a large r; n is not u nreasonable, ine pi 
iron Ka emission line near 6.4 keV (see, e.g., IPonti et al.l 120061 ; iMiller et al.ll2010l ). whose 
origin is generally attributed to hard X-rays irradiating the relatively colder accretion disc 
([Reynolds fc Nowakll2003l ; lMillerll2007l ). is quite narrow and does not show strong relativistic 
broadening, further suggesting that the accretion disc may be truncated at a large radius. 
However, by comparing the observed variability in the r eflection fea t ures s een in the X-ray 
data with reflection models that include light-bending, IPonti et al.l (120061 ) have suggested 
a much smaller inner edge radius (within a few r g ). We would like to point out that our 
estimation of the disc's r in assumes that the optical, UV, and FUV light are entirely coming 
from the accretion disc. However, ongoing star formation, as evidenced by the PAH lines in 
the Spitzer spectrum, may contribute some fraction (although most likely not a large fraction; 
see below) of the light that we are currently attributing to the disc. Any (presumably 
small) contribution from star formation will reduce the disc contribution, which in turn will 
likely imply a smaller disc-truncation radius (assuming that T in of the accretion disk, which 
determines the peak of the optical/UV SED, is not changing). 



An additional contributor to the optical and UV flux is reprocessing of the more energetic 
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EUV/X-rays from near the central black hole. A recen t work studying the long-term optical 



and X-ray light curves of NGC 4051 bv iBreedt et all (120101 ) hints at a lag of 1.2+^° days 



with the optical lagging behind t h e X-r ays, thus favoring a reprocessing origin. On the other 
hand, results by IShemmer et all (120031 ) suggest that p art of the optical lea ds the X-ray by 
~2.4 days. Based on their detailed reprocessing model IBreedt et all (120101 ) concluded that 
there seems to be no common origin for the observed X-ray /optical variability, and that not 
all of the observed optical variability can be attributed to reprocessing. Contribution by 
reprocessing would tend to lower the disk r in . Therefore the r in deduced here should perhaps 
more realistically be seen as an upper limit on the actual inner radius. 

An independent measurement of the inner radius may come from a self-consistent mod- 
eling of the soft excess in this source. While the true nature of the soft excess is still 
Galloir 



debated (see , e.g. 



2006 



Middleton et al.ll2007t iFabian et al.ll2009h . a survey of 34 type 



1 AGNs by ICrummy et al. J2006h a nd det ailed spectral and t iming analysis of the NLS1 
galaxy 1H0707-495 by IFabian et aP fj2009h and IZoghbi et all (boioh showed that the soft 



excess can be well modeled by relativistically blurred photoionized disc reflection model. 
Assuming a maximall y rotating Kerr black h ole, disc reflection model fits to XMM-Newton 
data of NGC 4051 by I Crummy et all (120061 ) give an inner disc radius of ~1.24r g . A more 
physical treatment of the scenario, where the X-ray emission from the jet-base (which can 
act as a source of irradiation) is self-consistently used to compute disc reflection features, 
would help in understanding soft excess seen in the source. If the inner disc radius is indeed 
close to the innermost stable circular orbit (ISCO) as suggested by the reflection model fits, 
this would require most of the flux, which is attributed to the accretion disc in this model, 
to have a non-disc origin. 



4.3. Stellar contribution in the optical/UV 



The HST/STIS spectrum given in figure 6 of ICollinge et all (120011 ) suggests that the 
continuum flux between 1200-1700A is approximately constant at F^~2xl0~ 14 erg cm -2 
s _1 A -1 , implying an observed isotropic luminosity L o b s =1.2xl0 41 (D/10 Mpc) 2 erg s^ 1 . In 
order to estimate the contribution from stellar flux at these wavelen gths we used the synthetic 
SED s of simp l e stel lar populations published by iMarastonl (120051 ). We used a model SED 
with lSalpeterl (119551 ) initial mass function (IMF), Solar metallicity, and red horizontal branch 
morphology (HBM). Since we are mainly interested in the blue/UV photons, the main factor 
that affects the luminosity is age, and influence of other factors (IMF, metallicity, HBM) 
are comparatively much smaller. For the chosen set of parameters, the IMarastonl (120051 ) 
models give SEDs as a function of the age of the stellar population. Given the shape of 
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the STIS continuum, a relatively young age of 



Maraston model between 1200-1700 A gives a model luminosity of L rno H e i=3.25xl0 33 - , - • 1 



^0.1 Gyr is needed. Integrating the 0.1 Gyr 
^modei^-^^J-u ergs" 

Therefore L b s /L moc iei ~ 3.7 x 10 7 . Since the Maraston models are normalized to 1 M , this 
implies that ~37 million young stars would be needed within the STIS aperture in order to 
have the same luminosity. Such a stellar concentration, while not implausibly large, is still 
at the hi gh end of what h as been observed in globular clusters and in nuclear star clusters 



(see, e.g. 



Cote et al.ll2006h . 



However there is a very important issue here that we have not considered: in reality it is 
well known that the stars in globular clusters are quite old, and so are the stars in the bulge. 
Nucle ar star clusters t ypically consist of a somewhat older stellar population. According 
to the iMarastonl (120051 ) models, the number of stars needed to produce the STIS observed 
luminosity would have to be increased by 4 orders of magnitude if the stellar population was 
1 Gyr old, inste ad of the as s umed 0.1 Gyr. On the other hand, recent results from Spitzer 
spectroscopy by ISani et al.l ( 120101 ) show strong PAH emission at 6.2/im from narrow line 
Seyferts, indicating ongoing star formation within the circumnuclear region included inside 
the IRS's slit-width (3.6"). Clearly, high-resolution MIR observations would be required to 
better understand ongoing star-formation in the nuclear region of NLSs. The variability 
seen in the optical wavelengths however argues against a scenario where stellar emission 
dominates the total observed optical flux. Based on our modeling, <10% of the Spitzer 
observed flux is due to the jet. The jet contribution is higher in the 0.1" Gemini spectrum, 
but still appears to be dominated by sources other than the accretion disk and the jet. 



The results of modeling the broadband SED of NGC 4051, along with similar results 
from modeling stellar mass black holes as well as low-luminosity AGNs like Sagittarius A* 
and M81*, broadly suggest that the underlying physical mechanism driving accretion inflows 
and outflows do not depend on the mass of the accretor and rather most of the observed 
properties scale with mass. Evidence for the presence of this mass-scaling is not only seen in 
broadband spectral modeling and the existence of the fundamental plane, but also in various 
other recent observations which include (but a re not limited to) detection of quasi-periodic 
oscillations in the NLS1 g alaxy RE J1034+ 396 dGierlihski et al.l hood: iMaitra fc Mfllerlboiol : 



Middleton &: Pond 120101 ) , variability ( see iMcHardvl l2010i for a recent review) , frequency- 



dependent lags between energy bands ( Arevalo Sz Uttley 20061 and references therein), iron 
line a nd other associated disc reflection features (see, e.g., iReynolds fc Nowak 2003 ; Miller 



20071 ). It is interesting to note that recent X-ray and optical observations of gravitationally 



lensed quasars RX J1131 1231 and HE 1104-1805 suggest th at the X-ray emitting region of 



these quasars are quite compact (~10r g or smaller; see, e.g., IChartas et al.l 12009 



Dai et al. 



20101 ). As our modeling of NGC 4051 presented here shows, such a compact region of emission 
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is expected if the X-rays are produced at the base of a jet. 

Typically the inner disc radius is close to that of the ISCO in XRBs in soft state, 
whereas the inner-disc radius we obtain here is much larger and more characteristic of hard 
state XRBs where the disc is truncated (although the radius we obtain is most likely an 
upper limit as discussed above, but otherwise it would require a rather improbably high 
number of extremely young stars to produce the observed UV flux). Based on our modeling, 
the continuum luminosity of the accretion disk is 6.1 x 10 41 erg s _1 (~0.3%LEdd) and the 
isotropic luminosity from the power law component between 10 16 -10 19 ' 5 Hz is 1.33 x 10 42 erg 
s _1 (0.6% LEdd> about twice the disk flux). Integrating the broadband model flux between 
10 9 -10 19,5 Hz and assuming isotropic emission gives a luminosity of 2.3 x 10 42 erg s _1 or ~1.1% 
L Edd , in good agreement with the bolometric luminosity obtained with the L ho \ & 9\L\ (5100 
A) relation. Therefore, while the timing properties of NGC 4051 are similar to soft state 
XRBs, the SED and luminosities are more akin to hard state XRBs. Thus it is unclear if 
the XRB state definitions can be straightforwardly carried over to AGNs, or if perhaps NGC 
4051 is in a transitional phase between hard and soft states. 

It would be interesting to test whether outflow-dominated models can explain the SED 
of other Seyfert galaxies (Maitra et al. in prep.). As shown in this work, isolating the 
AGN emission from starburst, dust etc. is one of the biggest challenges for this project; 
and it is only in the very nearby Seyferts that we currently have the required resolution and 
sensitivity to see these relatively weak jets. The enhanced sensitivity of the EVLA at radio 
frequencies, and ALMA in millimeter and submillimeter wavelengths will soon provide data 
with unprecedented sensitivity and coverage, zooming into the heart of the central engines of 
the nearby AGNs to study accretion/ejection processes in much greater detail than before. 
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Fig. 1 — Broadband SED of NGC 4051. Chandra/ A CIS (solid black lines) and VLA 8.4 
GHz (red points) data t aken from iKing et al.l (1201 it ) . Suzaku/PIN data (magenta points 
connected by lines) from iMiller et al.l ( 120101 ). Swift/BAT data (brown points conne c ted b y 
line) from publicly available NAS A archives. EUVE data (violet) from lUttley et al.l ( 120001 ). 
FUSE data (blue poi nts) are from lKaspi et al.l ( 2004 ). HST/STIS fluxe s (magenta line) from 
Collinge et al.l ( 200ll ). an d 5100A continuu m data point (orange) from Dennev et al. ( 2009 ). 
The Gemini/NIFS data jRiffel et al.l l2008h and the Spitzer /IRS data ISani et al.l feoioh are 
shown in cyan and green respectively. The SED was corrected for Galactic absorption 
assuming a column density of E(B-V)=0.013 and Nh=1.3x10 20 cm -2 in the direction of 



NGC 4051 (lElvis et al.lll989h . and assuming the ext inction law o f ICardelli et all (11989( 1 . The 
fits to high-resolution Chandra spectra discussed in IKing et al.l (120111 ) suggest that intrinsic 
absorption is negligible. The thick solid lines show the various model components as well as 
the total model predicted SED. 
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Fig. 2. — The solid black line in each of the above panels is the best-match model SED shown 
in Fig. [TJ The dashed blue lines show the SED when one of the jet model parameters (labeled 
near the top-right for each panel) was varied by ±10%. The remaining parameters were fixed 
to the best-match model. The dotted brown lines show the corresponding variability in the 
SEDs when the model parameter is changed by ±20%. 



Table 2. Comparison of jet model parameters for different sources. The values for NGC 
4051 are b ased on this work, an d the others taken from a compilation of prior published 
literature by lMarkoff et al.l ( 120111 ). An asterisk (*) next to a parameter value indicates that 
the parameter was not varied while searching for the best-match model. 



Parameter (unit) NGC 4051 HS-XRBs Sgr A* M81* 



Nj(L Edd ) 


1.7xl0~ 2 


10~ 7 - 10" 1 


io- 8 


10~ 5 


r o(r g ) 


4 


2-100 


2.5 


2.4 


ho/r 


1 


0.4-1.7 


0.4-2.5 


5-15 


Tc(K) 


4xl0 10 


2 - 5 x 10 10 


10 11 


10 11 


k 


1* 


1-5 


>10 


1.4 


Zacc(^*g) 


5 


10-400 


>10 4 


144 


V 


2.7 


2.4-2.9 


>3.8 


2.4 



